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Abstract 



We study the production of charmonium hybrids in B decays. We use the 
operator product expansion and nonrelativistic QCD to organize the various 
contributions to this process. We express the decay rate in terms of a few 
matrix elements which eventually will be fixed by experimental measurements 
or calculated on the lattice. While the Fock state expansion is problematic 
for hybrids, in that there is no perturbative large mass limit, the operator 
product expansion still provides a model independent framework for phe- 
nomenological calculations of hybrid production and decay. We then use a 
simple flux tube model to estimate the branching ratio B — > ip g + X, where 
ipg is a J FC = H hybrid, the large production of which could help resolve 
the low charm multiplicity and semileptonic branching ratio observed in B 
decays. We also investigate the possible contributions of tensor charmonium 
production. We observe that it is unlikely for either effect to be large enough 
to play a significant role in resolving these problems. 
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I. INTRODUCTION 



The hybrid states remain an enigmatic feature of hadronic physics. Their existence is 
predicted in QCD (at least in the large N limit jl|), yet most of their properties have been 
studied only within models or on the lattice. This is true even of quarkonium hybrids, 
because, unlike ordinary quarkonium, hybrids are characterized by excitations of the non- 
perturbative gluon field. Even for large mass tuq, there is no reason for a QQ hybrid to be 
a compact object, and neither the lattice nor models support such a picture. The absence 
of a limit in which hybrids are well understood makes a rigorous analysis of properties such 
as their masses, decays, and production rates difficult if not impossible to achieve. 

Yet it is still important to do what one can. There is experimental evidence that hybrids 
exist in the light hadron spectrum H. As for charmonium hybrids, it has been suggested 
that they might play an interesting role in nonleptonic B decays, in particular, in resolving 
the puzzle of the low semileptonic branching fraction . Essentially, the problem is that there 
is a discrepancy between the CLEO experimental measurement [[J 

ggf = ~* Xlu) = 0.1019 ± 0.0037 (1) 

1 tot 

and theoretical estimates which prefer > 0.115 [|]||. Here r tot = T(b — > cc + X) + T(b — > 
c+X)+T(b — > no charm) is the total decay rate of the B meson. Assuming that T(B — > Xlv) 
can be predicted accurately, the simplest resolution is to boost the nonleptonic decay rate, 
typically either by enhancing Tib — > cc+X) or T(b — > c+X) within the Standard Model [§],§, 
or by enhancing T(b — > no charm) with new physics M. The first option is constrained, 
however, by the measured charm multiplicity from CLEO ||, 

n c = 1 + nb^cc + X) _ T(b -> no charm) = ^ ± Q ^ ^ 
Ttot Ttot 

which apparently cannot accommodate a large rate T(b — > cc + X). While new physics 
remains a possibility, it is certainly more conservative to look for new mechanisms to resolve 
the problem within the Standard Model.B 

One such possibility is the suggestion that a significant fraction of the b — > cc + X 
transition is seen as charmless b — > s processes 0. If this were true, then the con- 
ventional estimate of BiB — > no open charm) = 0.052 ± 0.011 could be boosted to 
BiB — > no open charm) w 0.16. The mechanism proposed to account for this involves 
the production of a charmonium hybrid which, for dynamical reasons, does not decay to 
open charm. For example, in some models of hybrids there are "selection rules" which 
suppress decays into mesons with "similar" wavefunctions, or the hybrid can have exotic 
quantum numbers which prohibit its decay into certain hadronic final states. In this paper, 
we will investigate the production of hybrid charmonium in B meson decay. We will not 



Measurements from LEP of = 0.1112 ± 0.0020 and n c = 1.20 ± 0.07 are consistent with 
the Standard Model, so there are experimental issues to be resolved as well. 
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have anything new to say about whether the hybrids under consideration do, in fact, have 
the desired masses and branching ratios; this is an extremely model-dependent question, and 
one which ultimately must be answered by experiment. Rather, we will address the issue 
of whether, if such hybrids exist, their production rate could be large enough to make the 
proposed mechanism feasible. In doing so, we will identify the model independent aspects 
of the analysis and develop them as fully as possible. 

We will begin in Section II with a brief discussion of the properties of hybrids. Section 
III contains a general analysis of hybrid charmonium production in B decays, based on the 
operator product expansion and nonrelativistic QCD. In Section IV, we will use a simple 
flux tube model to estimate the matrix elements which govern the production of this state. 
For concreteness we will focus on charmonium with J PC = +_ , an exotic meson which 
cannot decay to D^*'D^*> . Sections V contains observations concerning tensor charmonium 
production, and we summarize in Section VI. 



II. QUARKONIUM HYBRIDS 

A hybrid meson is a resonance which contains, in addition to a valence quark and an- 
tiquark, excited gluonic degrees of freedom. Some hybrids have exotic quantum numbers, 
such as J PC = + ~ or which cannot be realized by a pure QQ configuration, and such 
states, if observed, are clearly hybrids (or four quark states) of some sort. Other hybrids 
have conventional quantum numbers, and hence can mix with ordinary mesons; here the 
distinction between the two is less clear. The gluonic excitation itself is realized variously 
in different models. For example, one may introduce into the constituent quark model a 
constituent gluon which couples to the QQ state to form a hybrid 0. Alternatively, bag 
models treat the motion of the diquark state in the background of a gluonic standing wave 



in the cavity [10]. Flux tube models realize gluonic degrees of freedom by introducing col- 



lective phonon excitations on the color flux tube which connects the quark and antiquark at 



distances of the order of A QCD |T]]JT2|JT3[[ . Hybrids also have been studied using QCD sum 



rules |14]] and on the lattice |15|JI^| . 



Despite the wide variation in these physical pictures, the models of hybrids display some 
common features. They typically predict the same quantum numbers for many of the lowest 
lying states, even if the masses are somewhat different. It also has been noted that the states 
in these models often possess symmetries which forbid or suppress hybrid decay into two 
mesons with identical wavefunctions. This particular property, if true, might prove useful 
for resolving the semileptonic branching ratio and charm multiplicity problems in B decays, 
since the production of charmonium hybrids which do not decay to £)(*)£)(*) could provide 
channels through which a ccs final state could hadronize and then appear as a final state 
with no open charm. Such a mechanism would allow for an enhancement in T(b — > cc + X), 
but with no enhancement, in fact a suppression, of n c . 

What would be required for this mechanism is one or more narrow hybrid charmonium 
states favoring light hadronic decay channels. At this point, the evidence for the existence 
of such a state is very model dependent, and has in no way yet been confirmed by genuine 
QCD calculations. For the purposes of this paper, we will posit that such a state exists 
and investigate its production in B decays. To be specific, we shall study the exotic hybrid 
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state with J PC = + ~, where strong decay to /}(*)/}(*) final states is forbidden by charge 
conjugation or parity. This state cannot be realized by a cc combination and therefore 
has a unique signature of a hybrid. However, the formalism developed here can be applied 
generally to the study of any hybrid charmonium production in B decays, as well as to 
hybrid photoproduction and hybrid quarkonium annihilation. 

A times, it will be useful for us to have a simple but physical model of hybrid mesons, 



and for this purpose we will turn to the original flux tube model of Isgur and Paton [pTT 
(This model has been elaborated considerably; see, for example, Ref. [0.) In this model, 
the quark and antiquark are separated by a distance of order 1/Aq CD and connected by a 
string of color flux. Transverse oscillations are then excited on the string, corresponding 
to collective degrees of freedom of the glue. The excitation state of the string governs the 
effective potential of the Q and Q. For example, with one phonon in the lowest excited 
state, the potential energy as a function of the separation r is 

V( r ) = -— + c + br + - [l - exp (-fb 1/2 r)] , (3) 
3r r 1 v /J 

where b ~ 0.18 GeV 2 is the string tension, c ~ —0.7 GeV 2 is a constant, and / is a parameter 
of order one which governs the r — > behavior of the flux tube. Already in this sector, a 
rich spectrum of hybrid charmonium states is predicted, with J PC = ±Zf , 1 ±=F , 2 ±=F and 
1 ±± . Of these, the +_ and 1 h states are exotic. The CP quantum numbers of these 
hybrids are related to the total spin S of the quark and antiquark by rjcTjp = (— l) s , as 
compared to rjcrjp = (— for ordinary mesons. Thus in this model the hybrids can mix 
with conventional states only via hyperfine transitions. 

Hybrids are different from ordinary quarkonium in an important respect which makes 
their study problematic, which is that they have no simple limit as mg — > oo. In this limit, 
ordinary quarkonium admits a controlled Fock state decomposition, in which the leading 
term is a QQ in a color singlet with fixed quantum numbers 2S+1 Lj. Higher order terms, 
such as color octet configurations, are suppressed by powers of a s (mq) ~ v, where v is 
the relative velocity of the Q and the Q. The expansion is possible because for large mq 
the quarkonium is a small state, of size (rnqv)^ 1 , whose interaction with the color field is 
governed by a multipole expansion. By contrast, the minimum of the potential @ is at a 
separation r of order 1/Aqcd, independent of mg. For large ttiq, fluctuations about r are 
small, but the state itself is not compact.! This appears to be a generic feature of hybrid 
models; for example, in a constituent gluon model, the Q and Q are in a color octet and 
repel each other at short distances. At large m,Q, one expects a situation somewhat like a 
heavy rigid rotor, with nearly degenerate rotational bands of states. This behavior has been 
observed in lattice studies |TB| . 



The absence of a controlled Fock state expansion for hybrid quarkonium can be under- 
stood from simple physical considerations. In the case of ordinary heavy quarkonium, a soft 



2 One also might imagine quarkonium hybrids with the structure of a small color singlet quarko- 
nium bound to a large glueball. There is little evidence at present for such bound states, and we 
will not consider them here. 
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glue configuration has a small overlap with the compact two-quark state, since the Compton 
wavelength of the soft gluons, of the order of A.qc D , is much larger than the distance between 
the quarks. By contrast, the size of the heavy hybrid state remains finite, of order Aq CD , 
in the heavy quark limit, thus allowing for a significant interaction with the soft gluonic 
modes. Yet, in this limit heavy quarks still become nonrelativistic, and thus the production 
of virtual heavy quark pairs is strongly suppressed. 

Nevertheless, there is still a sort of Fock state expansion for quarkonium hybrids. Al- 
though the gluonic degrees of freedom are intrinsically nonperturbative, we can organize the 
various components according to the quantum numbers of the Q and the Q. For example, for 
the +_ charmonium hybrid state which we will study later, we can write a decomposition 
of the form 

|W> = A\cc( 3 S 1 ) h8 + 9l ) + B\cc( 3 Pj) 1>8 + g 2 ) + C\cc( x S Q ) 1>B + g 3 ) + . . . , (4) 

where the subscripts 1,8 indicate the color of the cc pair and the Qi are various gluonic 
configurations. While one is always free to do this, here there is no model-independent 
hierarchy among the coefficients. This is related to the lack of a unique cc "baseline" 
state, and the lack of a multipole expansion to govern transitions between the different Fock 
components. (Note that heavy quark spin symmetry still governs hyperfine transitions, such 
as between cc( 3 S'i)i i 8 and cc( 1 Sq)i ) s.) Nonetheless, within particular models it is typically 
the case that a given hybrid is dominated by one or two cc configurations. For large m c , the 
c and the c move nonrelativistically in the potential V(r), and an expansion in powers of the 
velocity v is available. Hence it is useful to organize our calculation of hybrid production in 
B decays with a decomposition such as (^) in mind. 

III. HYBRID QUARKONIUM PRODUCTION IN B DECAY 

A. Factorization 

We turn now to the calculation of hybrid charmonium production in B decays, using the 
formalism of nonrelativistic QCD (NRQCD), an effective field theory useful for describing 
nonrelativistic degrees of freedom fll7| . The applicability of this technique rests on the fact 
that even though a hybrid is not necessarily small as tuq — > oo, the heavy quarks are 
still moving slowly and virtual QQ production is suppressed. In the usual way, we will first 
calculate the production rate in the framework of full QCD, separating hard and soft degrees 
of freedom. Then we will compute the rate in NRQCD, and use the perturbative calculation 
to match the coefficients of the NRQCD matrix elements. 

We begin by writing the inclusive B meson decay rate to a (hybrid) charmonium state, 
B — > ipg + X, in the B rest frame 

rQ0D = i / S n / jSm^ &- p -z *) i ^ + *i*»»' b >i 2 • < 5 > 

where P is the momentum of the ip g and hi are the momenta of the other particles in the 
final state. Note that we keep the ip g off mass shell for now; we will impose this condition 
later, when we match onto NRQCD. The effective Hamiltonian takes the form 
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n eff = ^v cb v; s [d^OiOi) + cmom] 

01 = ^(1 + 75)^' c j7 M (l + 7 5 )&\ 

02 = ^(1 + 75)^^(1 + 75)^', 



(6) 



where Ci(m&) ~ —0.29, ^(ni;,) ~ 1.1, and we have neglected small contributions from 
penguin operators. Note that the scale dependence of the operators is compensated by the 
scale dependence of the Wilson coefficients, which we will no longer make explicit. Using 
Fierz identities for the color matrices, Eq. (f|) may be rewritten to bring it to a form more 
suitable for our calculation, 



G F 



(Ci + |C 2 ) cr M (l + 75 )c s 7 "(l + 75 )6 
+ 2C 2 c 7m (1 + 7 5 )T a c sY(l + 7s)T a 6 



(7) 



with T a = A a /2. Here we are interested in the production of the hybrid charmonium states, 
in which there is no preference for the color singlet cc combination. Hence the first term in 
Eq. (|7|) may be considered as a perturbation of the second, since (Ci + |C 2 ) <C C 2 , and we 
will neglect it. Defining r M = 7/Lt (l + 75), we may write 



N-l 

n 

i=i 



d 3 L 



(2tt) 3 2E, 



(2tt) 4 J( 4 > 



2G^|K fe \/ c ;| 2 C , 2 2 ^ {B\bY,T a s\X l ){X l \sY u T b b\B) 



: (0|cr^T a c|^ + X 2 )(^ g + X 2 |cr"T 6 c|0) 



2G 2 P iKtKtl 2 O 2 £? 6 .C„T • 



fj,v ab 



(8) 



Using the optical theorem, and the fact that all asymptotic states are color neutral, B 
may be related to the imaginary part of a time-ordered product of color singlet currents 

S ab 



ab 



Bf v = -Im ? J d*ye* p y(B\T{jl(y),j u (0)}\B). 



(9) 



This expression for B^ v could be a starting point for the usual calculation of l/ra& corrections. 
The time ordered product (||) may be expanded in inverse powers of nib, yielding a tower 
of operators whose expectation values in the B meson are organized according to the heavy 



quark expansion ||18||. Here we will concern ourselves only with l/m c corrections and use 



the lowest order result, which coincides with the parton model expression for B 



ab 

)1V1 



B 



ab 



^rt(^-P) 2 ) 



(I \ 2 
- 9 to, 1 2~Pb ■ P M aJVPta/ 



HIT 



mi 



1 



{PpPbv + PuPbp) + -^ZpvapPn 



mi 



(10) 
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where we neglect the mass of the s quark. The leading nonperturbative corrections to this 
result would be proportional to Ai^/tti 2 ,. Neglecting these small effects, the decay rate in 
full QCD is then written as 



■pQCD 



N 
2rrih 



d 4 P 



(27T 



(11) 



with iV = 2Gp \V c bV* s \ 2 C\. The expression of Eq. ([LI]) must then be matched with the 
equivalent expression in NRQCD, 



■pNRQCD 



N 
2rrih 



d A P 



E 



rrv 



dn-4 



(0^)27r5(P 2 - 



m 



(12) 



where 0% 3 represents a set of NRQCD operators of increasing dimension d n , and we have 
extracted the on-shell 5 function explicitly. We are unable to calculate matrix elements of 
these operators in a model independent fashion, so as usual they are left as free parameters. 
Note that the same matrix elements would govern the total annihilation width of the hybrid 
or the hybrid photoproduction cross section. Eventually, one might hope to extract the 
leading matrix elements from a small number of experiments. 



B. Perturbative QCD Calculation 

We now turn to the perturbative calculation of , 

C% = (0\cT»T a c\^ g + X 2 )(ij g + X 2 \cT»T b c\0) = (Q\cT»T a cP- cc cT»T b c\0), (13) 

where we have introduced the projection operator P 5c = \cc + X 3 ) (cc + X 3 \ onto the desired 
cc state within the ij) g . The relativistic four- component spinors u c may be boosted to the cc 
center of mass frame and then replaced by nonrelativistic two-component spinors £ and r), 
using the identities! Jl7j 



u c ^T a u c = 



2E£a j T a i 1 - 



E c + m. 



u cl ^T a u c = --2 P^T a v + 2iA?£+(q x 



E, 



(14) 



where is the Lorentz boost matrix. For clarity, we introduce primes for one of the pairs 
of spinors. Keeping terms up to order q 2 , we obtain at lowest order in QCD perturbation 
theory, 



C a T = 4m 2 A^< 



1 + 



q 



12 



2m? 2m? 



i l] a j T a r]'rfa k T h i 



3 Our definition of 75 differs by a sign from that of Ref. [17|. 
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1 



2m 2 c 
1 

m] 



+ -C' t (q'xa)^V(qx ff )¥e 



(15) 



+ P M P v 



q' 2 + q 2 



2m;? 2ml 



where there is an implicit projection onto the desired spin and angular momentum configura- 
tion of the quarks. Note that Eq. (|I5D , along with (Q) and ([11]), is valid for any charmonium 
hybrid state. If we now specify to the case of unpolarized hybrids, we arrive at a simple 
expression for r^ CD , 

r QrL) = JL I d ! P . 27rml (1 - + 2^)5((p b - Pf) x Ami 



2m b J (2vr) 4 



+ 



./2 



N - ( 1 - - + o q I ^T a ri'^T a i 



(16) 



2m? 



f'V(q' • a)T a ri'rfo j T a ^ + f V'TV^V (q • a)T a £ 



where /i c = m 2 c jm\ and /i^, = P 2 /ml, in anticipation of matching onto the hybrid state for 
which P 2 = m\. Eq. ([H]) expresses the decay rate in terms of matrix elements of operators 
of increasing dimension, and the projection onto an unpolarized hybrid is understood. The 
explicit factor of Am 2 comes from the original normalization of the relativistic cc state. 

C. NRQCD Calculation 



To perform the matching, we have to calculate T(B — > ip g + X) in perturbative NRQCD. 
The expansion takes the form 



■pNRQCD 



N 



d 3 P 



2m b J (2tt) 3 2^ 



2m„ 



cl(" scopes,)) + cics )(o^cs )) 



+ 
+ 



m 



mt 



mt 



mt 



where the NRQCD operators are given by 

(OpfSx)) = (0|xVT>P^VV'T a X |0), 
(OPCSo)) = (0\ X j T a i;P^T a x\0), 



1 



(Ot 3 ( 3 Pi)} = -(0|x t (-fDxayT^^t(_i D x ( T)^ X |0), 



(17) 



/ r>i>g 


v = 


i 

2 




/" T~\\2 r pa,j, z~> n fA ^.irncb* An\ i /n U ,t *i r Pa„i, ID n lA ~i if i~f~\ 


) 1 X|0) 


(pi- 


^o)> = 


1 
2 


;<oix f (- 


f D) 2 T>P^t T « x |0> + (Olx^VP^H D) 2 T a X |0>~ 




(Qi 3 


: 3 ^i)) = 


1 
2 


;(oix f (- 


f D)'(-| D • a)T>P^Vr a X |0> 





-D-cr^xlO) 



(18) 



Here P^, 9 is a projection onto the hybrid state ip g , and D = D — D. By convention, in Eq. ([] 
the state \ip g ) is normalized nonrelativistically. This is made explicit by the prefactor 2m^, 
so the coefficients Cj remain dimensionless, as they are in Eq. (0). Color singlet operators 
could also be included. Note that each derivative insertion brings a power of the velocity of 
the heavy quark, p/m c ~ v. For instance, Og 9 ( 3 Pi) is suppressed compared to 0$ 9 ( 3 Si) by 
v 2 . Nonrelativistic kinematics clearly is crucial for controlling the number of nonperturbative 
parameters, and it is therefore not practical to extend these calculations to the case of lighter 
hybrids. 

It is now straightforward to evaluate the NRQCD coefficients, by matching full NRQCD 
states to perturbative ones. We have kept the normalizations of the states explicit so that 
we can treat them carefully. With our conventions, we must match quantities in the manner 

lml£*T a rtrfT a t <— > 2m^ (0\ X j T a ^P^T a X \0) ■ 2n5(P 2 - m%) , 

where both sides have mass dimension two. We find 



(19) 



°8 



Si) 



ClC^o) 
Cl( 3p i) 

cf( 3 sl) 
cteso) 

CK^) 



2-nm\ 



SfJ-ip C 1 ( 3 S ) 
4// c (l + 2/v) 8 
2C 8 1 ( 3 S 1 ), 



^)(l + 2^)5(( Pb - Pf 



Si, 



-clC 

~ U 8 



50) - 

51) , 



(20) 



As discussed in Section II, the Fock space decomposition of hybrid states is not controlled by 
a multipole expansion. However, the number of matrix elements important for production 
processes still can be reduced by keeping those that are leading in powers of l/m c , or 
equivalently, in powers of v. For instance, in this case only two matrix elements appear at 
leading order. (We remind the reader that in the case of hybrids, the identification of the 
leading matrix contributions is unavoidably model dependent.) 

We can now perform the remaining integrations to obtain the partial decay rate for 



r. 



4(1-^)2(1 + 2^) 



(op( 3 s 1 )) + 



(OpCSo)) + — 2 (0^( d P 1 )) 

mi 



+ —(P+'?S 1 ))- 



4/i c (l + 2^) 

iTD^Bfl 



4/i c (l + 2fx^)m] 



-<^ 9 (%)> 



mi 



■(QPCSi)) 



(21) 
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where we have defined T = Cl\V c bV* s \ 2 G 2 F m\/lAA'K 1 which is the free b quark decay rate 
induced by the second term in the Hamiltonian (|7|), for m c = 0. Note that as rnj, — > oo for 
m c fixed, the branching fraction for B — > ip g + X falls as l/mf. The terms proportional 
to (0 8 ( 3 5'i)) and (0 8 ( 3 Pi)) agree with the corresponding terms in Ref. fll9l , in which the 
production of ordinary P wave charmonium in B decays is analyzed. The formula ( |2T|) 
determines the semi-inclusive decay rate of B into a given hybrid charmonium state in 
terms of parameters which must be fixed by other experiments. We emphasize that these 
parameters are universal] for example, the photoproduction cross section will be expressed 
in terms of the same matrix elements, which depend on the static, equilibrium properties of 
the hybrids. 



IV. MATRIX ELEMENTS IN A FLUX TUBE MODEL 

Now that we have a general formalism for the production of hybrid charmonium in B 
decays, we would like to return to the question of whether such a mechanism can account 
for the "missing charm puzzle" and low semileptonic branching ratio of the B. For this 
purpose, we will focus specifically on the J PC = H hybrid, an exotic meson which cannot 
decay to D^*'D^*\ Although dynamics may forbid or suppress the DD decay channel of 
some conventional hybrids ||13|| , the effect is model dependent, and we will restrict ourselves 
to this more robust case in which the decays are prohibited by spacetime symmetries. Of 
course, it is quite possible for the + ~ hybrid to be heavy enough that it decays to multibody 
states with open charm, since its mass is not known model independently. However, we will 
assume for the sake of argument that it is not so massive. 

We will estimate the matrix elements (|18D in the flux tube model of Isgur and Paton |TT| . 
We are forced to use phenomenological models because of the lack of a first principles 
understanding of the structure of hybrid mesons. Eventually, these matrix elements should 
be computed on the lattice. Our model dependent calculation serves two purposes. First, it 
allows us to use the expansion fl2T|) to obtain a rough estimate of the magnitude of hybrid 
charmonium production in B decays. Second, we can explore the extent to which, within 
the flux tube model, the expansion is well behaved. We will use the flux tube model in 
its least elaborate form, and with some additional assumptions, on the grounds that this 
approach is sufficient for addressing these relatively simple questions. 

As discussed in Section II, the flux tube model is a model of nonrelativistic quarks moving 
in the potential (|3|), which incorporates the effect of a phonon excitation on the flux tube 
connecting the quarks. Because angular momentum can be exchanged between the flux 
tube and the quarks, the ground state of the Schrodinger equation is not an eigenstate of 
the quark quantum numbers 2S+1 Lj. Rather, for the H hybrid, it is a linear combination 
of 3 Si, 3 Pi and 3 Di states. In this model, 3 Si state has the lowest energy, with the 3 Pi state 



excited by approximately 0.24 GeV. The analysis of Isgur and Paton |]TT| finds a hybrid 
mass roughly in the middle, indicating significant mixing between the two states. The 3 Pi 
configuration is excited by approximately another 0.2 GeV over the 3 Pi state; therefore its 
contribution to the ground state +_ hybrid is likely to be small, and we shall neglect it 
here. There is no l So component in the one phonon sector, because in this model such states 



have the wrong charge under CP |T1]. The decomposition then takes the form 
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|0 + -> = A\cc( 3 St) 1>s + gi ) + S|cc( 3 Pi) ll8 + g 2 ) . (22) 

In the case of ordinary quarkonium, it is possible to make a model independent identification 
of the leading (color singlet) Fock component, with the admixture of other components 
controlled by an expansion in powers of v. For hybrids, the identification of a "leading" 
component is model dependent. Here we will choose to truncate the Fock expansion at 
these two quark configurations, and, for reasons discussed in Section IIIA, we will take the 
quarks to be in a relative color octet. However, the constants A and B will remain free 
parameters, subject to the constraint |A| 2 + |i?| 2 = l. 

To calculate the desired matrix elements in this model, we need wavefunctions which are 
regular as r — * 0. Hence we will solve the Schrodinger equation separately for the states 
of definite quark orbital angular momentum L, and use these solutions to estimate matrix 
elements such as (0 8 ( 3 Si)) an d (0 8 ( 3 Pi)). We neglect the direct contributions to the matrix 
elements from mixing between the eigenstates of L, which amounts to working only to first 
order in the phonon term in V(r). Within this model we can account explicitly for the 
coupling of the quarks to the short distance process, but not for the coupling of the flux 
tube to the soft gluon background. We will parameterize by £, a nonperturbative parameter 
presumably of order one, the amplitude for the color octet current to excite the one phonon 
excitation on the flux tube. We remind the reader that, within the present context, we have 
little insight into the actual magnitude of £, so much larger or smaller values than £ ~ 1 are 
also possible. 

The flux tube, which itself carries color charge, has the effect of decorrelating the color 
states of the quarks at either end. We will take a statistical approach, noting that an 
uncorrelated quark-antiquark pair has a probability | to be found in a relative color octet. 
In addition, we have to project the total (spin and orbital) J = 1 carried by the quarks and 
the unit angular momentum of the flux tube onto a hybrid of total spin zero. This projection 
leads to another factor of | in the matrix elements. Including these three factors, which 
combine to ^-£ 2 , the matrix elements may now be written in terms of the nonrelativistic 
radial wave functions Rs(r) and Rp(r). The result is |T7]] 

(0 8 CS 1 )) = ^e\A\ 2 x ^\R S (0)\ 2 

(Os( 3 Pi)) = l 1 e 



\A\ 2 x 


2tt ' 
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\B\ 2 x 






2^ ' 



(P s ( 3 S l )) = ~^e\A\ 2 x - Re R* s (0) V 2 i? 5 (0) 



\A\ 2 x 


3 


Re 




2^ 




\A\ 2 x 


1 


Re 




2^ 





WSi)) = --e\A\ 2 x-Re [R* s (0) V 2 R s (0) . (23) 

01 Z7T L J 

There is a linear divergence in V 2 i?^(r) as r — > 0. However, we may use the Schrodinger 
equation to rewrite (P 8 ( 3 5i)) as 

WSJ) = ^ 2 \A\ 2 x A 2m [V(0) - E] \R S (0)\ 2 

Ol Z7T 

= 2m[V(0)-E] (OsCSJ). (24) 

Here it is clear that this divergence is proportional to (a s /r)(0 8 ( 3 Si)), so it is compensated 
by a linear divergence in the first radiative correction to the coefficient of (0 8 ( 3 Si)). Since 
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we have not included these radiative corrections in our operator product expansion, we will 
subtract this divergence explicitly [17 1. There is also a logarithmic divergence which we have 
subtracted implicitly, and which changes the matrix element by an amount of order one. 
Since our primary purpose in computing (P 8 ( 3 Si)) is to show that it is small, this ambiguity 
is not a problem. Solving the Schrodinger equation numerically, with potential parameters 
b = 0.18 GeV 2 , c = —0.7 GeV 2 and / = 1, and a constituent charm quark mass of 1.77 GeV, 
we find 



<0 8 ( 3 Si)> 



27tt 
4 



e\A\ 



(Os( 6 Pi)) = ^-e\B\ 

4 



0.58 GeV 3 
0.03 GeV 5 



27vr 
4 



W s Si)> 
(QsC Sl )) 

Since in this model (O 8 ( 1 5 )) = (Pg^So)) 

/i^) 2 (l + 



\A\ 2 


0.14 GeV 5 " 


\A\ 2 


0.14 GeV 5 " 



0, the expansion (|21|) becomes 



To 



16 1 



27nm 2 



r 0.58|A|^ + 0.101-Bp + 0.08|A|^ - 0.03|A|^ + . . . 



(25) 



(26) 



where we have taken m c = 1.35 GeV. The second, third and fourth terms are formally 
suppressed by v 2 compared to the first, and we see that within this model they are indeed 
relatively small. Taking = 4 GeV and m;, = 4.8 GeV we find 



r( jB ^^ + x)~(2xio- 4 )e 2 r c 



(27) 



The inclusive rate T(b — > ccs) is suppressed relative to r by a phase space factor of approx- 
imately 0.25. Taking A = 1, to be conservative, and the amplitude for flux tube creation 
£ ^ 1, we obtain 



T(B ^i/j g + X)<(lx 10~ 3 ) r(6 -> ccs) . 



(28) 



According to this estimate, then, the production of this hybrid charmonium is more than 
a factor of ten lower than that of typical ordinary charmonium states. Of course, there is 
considerable model dependence in this calculation, and we have considered the production 
only of a single hybrid. Nonetheless, given the tiny overall magnitude of the result, we 
conclude that it is quite unlikely that this mechanism could play an important role in 
"hiding" the charm produced in b — > ccs transitions. The addition of additional hybrid 
channels will not change this qualitative conclusion. 



V. TENSOR CHARMONIUM PRODUCTION IN B-DECAYS 

The decay into DD pairs is prohibited by parity not only for some exotic hybrids, but 
also for the ordinary D wave charmonium states if) {2 ) and ip(2~ + ) (respectively, the cc in 
a 3 D 2 and a l D 2 configuration). Perhaps the production of these tensor mesons, and their 
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subsequent decay to light hadrons, could also "hide" the charm produced in underlying 
b — > ccs transitions. The decay rates for B — > ip(2~ ± ) + X have been computed using 
the same NRQCD techniques developed for conventional quarkonium in Ref. ||20|| , with the 
results 



m - ccCDj) + x) = r„ + 



r(B -> ccC 1 ^) + x) = r 6(1 2 (o^ 2 (%)) , (29) 

where the approximation = 2m c was made. Note that coefficients of the operators 
(Os( 3 5'i)) and (Og( 1 5'o)) are the same as in Eq. (|2l|) . Using velocity scaling rules and 
data from Z° decays, the corresponding branching fractions were estimated to be B(B — > 
^(2~) + X) ~ 1.1 5.2 %, and B(B -> ^(2"+) + X) ~ 0.7 %. 

However, for this mechanism to be responsible for boosting T(6 — ► no charm), the 
branching ratios of ip(2 ) and ip(2~ + ) to light hadrons must be sufficiently large. Un- 
fortunately, this does not appear to be the case here. For the ip(2 ), one estimates 
B(V>(2— ) -> Xci + 7) ^ 0.64, B(^(2— ) -> Xc2 + 7) ^ 0.15, and B(V>(2~) -> 
J/i/jtcti) ~ 0.12 pl| . Assuming that the rest of the decays are into light hadrons, we 
obtain B(ip(2 ) — > light hadrons) ~ 0.09 and the combined branching fraction B(B — > 
V>(2~~) + X) x fi(^(2~) -> light hadrons) ~ (1 5) x 10~ 3 . The addition of the V(2~+) 
does not change the situation: taking B(ip(2~ + ) — ► /i c ( 1 Pi) + 7) ~ 0.8 as the dominant 
decay mode and assigning the rest to £>(■?/> (2~ + ) — > light hadrons) results in the combined 
branching fraction B(B -> V(2 _+ ) + X) x B(V»(2 _+ ) -> light hadrons) ~ (1 2) x 10~ 3 . It 
is clear that unless the estimates of NRQCD matrix elements are extremely inaccurate, this 
mechanism cannot play an important role in enhancing the charmless branching fraction in 
B decays. 



VI. CONCLUSIONS 

In the spirit of nonrelativistic quantum mechanics, we have developed a QCD-based 
expansion for the production of hybrid charmonium in B meson decays. Although no con- 
trolled Fock state expansion exists for hybrids, even in the limit of large m c , our treatment 
provides a useful phenomenological starting point for model calculations of hybrid matrix 
elements. In particular, it allows us to organize the various contributions to the production 
process in powers of the charm quark velocity v, since in most models the quarks are, in 
fact, nonrelativistic for large m c . Although we have treated explicitly only production in 
B decays, the same matrix elements which appear here will also govern other production 
processes, as well as decays of hybrid states. 

Based on the nonrelativistic expansion, we have used the flux tube model of Isgur and 
Paton to estimate the rate for the production of the exotic J PC = + ~ charmonium hybrid 
in B decays, and found that it is considerably smaller than the production of ordinary 
charmonium. Hence it is quite unlikely that this mechanism can play a significant role 
in resolving the difficulties with the semileptonic branching ratio and charm multiplicity 
observed in some experiments. We stress that we have employed a very simple model, with 
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the goal of obtaining an order-of-magnitude estimate only. As models of hybrids become 
more refined, and as lattice studies develop further, more reliable calculations based on this 
expansion will become possible. 

Finally, we have restricted ourselves to an estimate of the production of the + ~ exotic 
hybrid, because its decays to £)(*)£)(*) are prohibited by spacetime symmetries and it is 
likely to be among the lightest charmonium hybrid states. Nevertheless, its mass is not yet 
known, and it is certainly possible that other hybrids, perhaps even ones with conventional 
quantum numbers, may be more important in B decays. This is a question which could 
be addressed with a more comprehensive and up-to-date survey of available models, but we 
would prefer to wait until the spectrum of hybrid charmonium, and the decays of the various 
states, have been revealed by experiment. 
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